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Abstract 



We have measured the first and second moments of the hadronic mass-squared 
distribution in B ^ XJu, for Piepton > 1-5 GeV/c. We find {Mj. - M|,) = 
0.251 ±0.066 GeV^, ((M|. - (M|.))2) = 0.576 ±0.170 GeV^, where Md is the 
spin-averaged D meson mass. From that first moment and the first moment 
of the photon energy spectrum in 6 — > S7, we find the HQET parameter Ai 
(MS, to order 1/M| and Poal) to be -0.24 ± 0.11 GeV^. Using these first 
moments and the B semileptonic width, and assuming parton-hadron duality, 
we obtain \Vcb\ = 0.0404 ± 0.0013. 
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The heavy quark hmit of QCD |T[] is potentially a very useful tool for relating mea- 
sured inclusive properties in B meson decay, such as semileptonic branching fractions, to 
fundamental CKM parameters like Vcb and Vub- The expressions for inclusive observables 
are expansions in inverse powers of the B meson mass Mb [|]^- At order 1/Mb, the non- 
perturbative parameter A enters, and at order 1/M^ two more parameters, Ai and A2 appear. 
Intuitively, these parameters may be thought of as the energy of the light quark and gluon 
degrees of freedom (A), the average momentum-squared of the b quark (-Ai), and the energy 
of the hyperfine interaction of the spin of the b quark with the light degrees of freedom 
(A2/MB). The parameter A2 can be extracted directly from the B* — B mass splitting [Q. 
The other two parameters can be obtained from inclusive measurement or calculated theo- 
retically with techniques capable of handling non-perturbative effects, such as lattice QCD 

i- 

There are two problems associated with the interpretation of measured inclusive prop- 
erties, one associated with the convergence of the expansion, and another with the validity 
of the assumptions underlying the expansion. The inclusive observables are expansions in 
powers of 1/Mb, and at each order more non-perturbative parameters appear. By order 
there are three parameters and at order 1/M^ another six parameters, pi, p2, Ti - T4. 
Without good estimates for the additional parameters we must rely on the rapid convergence 
of the expansion. The other problem is the validity of the assumption of parton-hadron du- 
ality implicit in this approach, and its potential for introducing additional uncertainties not 
included in the present estimates 0. Thus, the experimental determination of A and Ai 
with several different methods is necessary to support the validity of parton-hadron duality 

i- 

Much interest has been raised by the possibility of estimating A and Ai using hadronic 
spectral moments in semileptonic B decays @-@]. In this Letter we report a measurement 
of the first and second moments of the distribution in the hadronic mass-squared in the 
inclusive semileptonic decay b civ. For this analysis, the leptons are restricted to the 
kinematical region Pi > 1.5 GeV/c. In particular, we report measurements of {Mx — M^)) 
and ((M|. -M^f), where M| is the mass-squared of the charmed hadronic system Xc, and 
Md is the spin-averaged D meson mass, 0.25Md + O.TSM/)* = 1.975 GeV. The theoretical 
expansion for these two observables has been carried out to order 1/M^ and order jSoal 
in the MS renormalization scheme [0,0. (Here /3o = (33 — 2nf)/3 = 25/3 is the one-loop 
QCD beta function.) We also report the second moment taken about the first moment 
rather than about M|,, i.e., {{Mx — (M|^))^), the mean square width of the mass-squared 
distribution. (The theoretical expansion for this is readily obtained from those for (M|^— M|)) 
and {{Mx — M|,)^).) We use the first moment, along with the first moment of the photon 
energy spectrum in 6 37 0, to obtain Ai and an improved extraction of Vcb from the B 
meson semileptonic width. 

The data used in this analysis were taken with the CLEO detector at the Cornell 
Electron Storage Ring (CESR), and consist of 3.2 fb^^ at the T{AS) resonance and 1.6 fb~^ 
at a center-of-mass energy 60 MeV below the resonance. The sample contains 3.4 million 
BB pairs. We select events containing a lepton - /i or e - with momentum between 1.5 
and 2.5 GeV/c. We "reconstruct" the neutrino in the event using energy and momentum 
conservation of the entire event, exploiting the hermiticity of the CLEO detector. The 
neutrino energy is taken as the difference of twice the beam energy and the sum of the 
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energies of all detected particles, while the neutrino momentum is the negative of the vector 
sum of the momenta of all detected particles. Considerable effort was expended to remove 
double counting between calorimeter and tracking chamber measurements. To assure a well- 
measured neutrino, we require: a neutrino mass consistent with zero; no additional leptons 
in the event; a measured net charge of zero for the event. The "neutrino reconstruction" 
aspect of this analysis is similar to that of Ref. [|10|, and is described in detail in Ref. [pT] . 



Event shape requirements are applied to distinguish the jetty event environment typical of 
e'^e~ — > qq light quark pair production from the more isotropic environment of e'^e~ BB 
events. We achieve a sample consisting of 89% e+e~ — > BB and 11% from the continuum, 
with an efficiency for the desired events of ~2%. The desired semileptonic B decays, h —>■ ciu, 
represent 95% of the e^e^ —>■ BB sample while the remaining consists of (2.8 ± 0.6)% 
secondary lepton production (from Monte Carlo simulation) and (2.1 ± 1.1)% b uiu 
(using \Vub/Vcb\ =0.07 ±0.02). 

We determine the mass of the hadronic system X in B Xclv from the lepton and 
neutrino momentum vectors alone: 

Ml = {Eb -Ee- E,f - {Pb -Pi- Puf 

= Ml + Ml - 2EBEe, + 2\PB\\Pi,\ cos0,,,b . (1) 

For B mesons produced at the T{4S), Eb and \Pb\ are known and constant, but the angle 
between the B and iu system varies from event to event, and is not known. Since \Pb\ is 
small (300 MeV/c), we approximate Mx by dropping the final term in Eq. |I|, writing 

m| = M| + Ml - 2EBEi,. (2) 

The background-subtracted M^ distribution, consisting of 11900 B meson decays, is 
shown in Fig. [l|. The background from continuum events has been subtracted using the data 
collected below the T(4S') resonance, scaled to the luminosity of the on-resonance data and 
corrected for the dependence of the production cross section on beam energy. The small 
backgrounds from secondary lepton sources and from b ulv decays, which we obtain from 
Monte Carlo simulation, have also been subtracted. 

For the purpose of extracting the moments of the M^ distribution, we divide the b c^v 
decays into three components: B D£u, B —>■ D*ii', and B Xh^i^ where Xh represents 
all the high-mass charmed meson resonances as well as the charmed non-resonant decays. The 
individual components are shown in Fig. |l|. We use measured form factors [0 to model the 
B Div and B D*£v decays. The true M|- distributions for B Diu and B D*iu 
are narrow resonances at M|, and M|,* . The widths of the Monte Carlo predictions in Fig. |1] 
for these resonances are dominated by neutrino energy-momentum resolution and our neglect 
of the last term in Eq. |T|. The high-mass contribution, B XhC-V-, is modeled using six 



resonances above the D* with the decay properties specified by ISGW2 form factors [13 



and also non-resonant multi-body final states such as B Dniu and B D*Tcii', which 



are decayed according to the prescription of Goity and Roberts |14 . 

A fit of the Monte Carlo to the data M^ distribution determines the relative contributions 
from B Diu, B D*lv and B X^lv. The relative rates and the generated masses are 
used to calculate (M|^ — M|,) and {{Mjr — M|,)^) of the true M\ distribution. Equation || 
shows the derivation of the average mass squared, Mx, from the relative rates. 



5 



2000 





— Fit 






1500 

CM 

> 

0) 

O 


D*£v 




If) 




O1000 




(0 

c 




> 

HI 

500 






IVI^(GeV^) 



10 



FIG. 1. Measured distributions, for background-corrected data (points), Monte Carlo 

(solid line), and the three components of the Monte Carlo: B — s- Dlu (dashed), B D*i!.v 

(dotted), B Xh^v (shaded). The normalization of each component is derived from a fit to the 
data. 



{Ml) = ro-Ml + rz). ■ M|,. + rx, ■ (M^^), 



(3) 



where r/j is the rate of -B ^ Dlv production compared to the combined rate oi B ^ Diu, 
B D*ih', B —>■ Xh^i^, and similarly for r^. and rx^- The individual values obtained for 
^D, fD* and rxff, while perfectly consistent with world average branching fractions [|15|, are 
not well determined and are sensitive to the model chosen for B —>■ Xh^v- The moments, 
however, are well-determined and stable against model changes, as discussed below. We find 
{Ml - Ml) = Ml = 0.251 ± 0.023 ± 0.062 GeV^ {{Ml - Ml)'^) = M2 = 0.639 ± 0.056 ± 
0.178 GeV^ and {{Ml - {Ml)f) = M2' = 0.576 ± 0.048 ± 0.163 GeV^ where the errors 
are statistical and systematic, in that order. The experimental errors on {{Ml — {Ml))"^) 
are somewhat smaller than for {{Ml — Ml)"^) and have a smaller correlation with the first 
moment. (A correction for final state radiation, not included in the Monte Carlo samples 
used in our fits, has been applied, using PHOTOS [T^].) 

The errors on both first and second moments are dominated by systematic errors. The 
leading contribution is from the simulation parameters that impact neutrino resolution: 
photon identification efficiency, tracking efficiency, and the rate of additional neutrals such 
as and additional neutrinos; it amounts to ±0.058 GeV^, ±0.140 GeV^, and ±0.129 
GeV^, for Ml, M2, and M2', respectively. 

The second leading source of systematic error is from the models for the high-mass 
contribution to the Ml distribution. We have varied aspects of the high-mass component in 
order to quantify the sensitivity. The six contributing mass states of the resonant component 
(above D*) have been systematically dropped singly, in pairs and in triplets so as to vary the 
internal structure of the resonant model. Taking the r.m.s. deviations of these variations, 
we find errors of ±0.015 GeV^, ±0.090 GeV^ and ±0.083 GeV^ for Ml, M2, M2'. 
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Another contributing uncertainty arises from the lack of knowledge on the amount and 
shape of non-resonant contribution to the high- mass component. Although we fix the frac- 
tion of non-resonant to resonant high-mass states during a fit, we systematically vary this 
fraction over the limits that the data allow. A one unit variation of fit determines a 
systematic variation of 0.011 GeV^ for Ml and 0.060 GeV^ and 0.054 GeV^ for M2 and 
M2'. Systematic errors other than those from neutrino resolution simulation, modelling of 
high-mass resonances, and modelling of non-resonant high-mass decays, such as the subtrac- 
tions for secondary leptons and for h — u^v, the final state radiation correction, and the 
B — > D^*^iu form factor uncertainties, are negligible by comparison. 

As an alternative to the default Goity-Roberts parameterization, we have also used a 



phase space model to generate the four-body non- resonant decays [|TT|. This phase space 
model generates, on average, higher mass states than the Goity-Roberts parameterization 
but yields hadronic mass moments consistent with those obtained from the Goity-Roberts 
parameterization. This observation emphasizes the fact that the data essentially constrain 
the product of the average mass squared and production rate while these quantities may 
individually vary significantly. 

The correlation coefficients between errors of first and second moments are positive, and 
substantial. They are +0.71 for Ml - M2 {+0.56 for Ml - M2') for the statistical error, 
+0.50 (+0.34) for the systematic error, and +0.52 (+0.36) for the total error. 

The expressions p|JT7[] for the hadronic mass moments in B ^ X^lv, to order (3oa1 and 
1/M^, subject to the restriction > 1.5 GeV/ c, are given in Eqs. ^ and ||. (Due to technical 
difficulties, the coefficients of the jj^^ terms were computed without the 1.5 GeV lepton 
energy restriction, and so are only approximate, beheved good to ± 50%.) 



^ ^ ^ ^' = [0.0272^ + 0.058/^0^ + 0.207#(1 + 0.43^^^) + 0.193^ + 1.38^^ + 0.203 „ 



2 



J^2 • ^ ■ r-.^. - Mgy -TT / ■ - —Ml, — —Ml, 



+ 0.19^ + 3.2% + 1.4|^ 

nig IVlg Ulg 

+ 4.3-£^ - 0.56-^ + 2.0^ + 1.8-i- + 1.7-^ + 0.91-^ + 0(1/M^)], 



^^^^ ^^hilll = [0.00148^ + 0.0025/3o4 + 0.027^^ + 0.0107^ - 0.12-,^ 

L 71- ' '^vr^ Mg TT Mg Mi 



+ 0.02^-0.06||-0.129|| 



B 



- 1.2^ + 0.0032-^ - 0.12^ - 0.36^ + 0(1/M^)], (5) 



In these expressions, Mb represents the spin-averaged B meson mass, 5.313 GeV. 

The 1/M^ parameters pj, Ti are estimated |0, from dimensional considerations, to be 
~ (0.5GeV)'^. Taking values of p2 and Ti through 7^ to be 0.0±(0.5GeV)^, taking pi (believed 
to be positive) to be |(0.5GeV)3±|(0.5GeV)3, taking A2 = 0.128 ± 0.010 GeV^ (appropriate 
with a calculation to order 1/M|) 0, and using ^^(mfc) = 0.220, the expressions combined 
with our measurements define bands in A — Ai space. The band for the first moment is 
shown in Fig. ^ The dark grey region indicates the error band from the measurement; the 
light grey extension includes the error from the theoretical expression, in particular from the 
pi — T4 terms and from the scale uncertainty (as(mfe/2) = 0.275 to as{2mi,) = 0.176). 
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In the preceeding Letter , we presented measurements of the first and second moments 
of the photon energy spectrum in 6 — 57, and gave the OPE expansion expressions for those 
moments, again vahd to order /3o«s and 1/M^. Again, equation plus measurement defines 
a band in A — Ai space. The band for the first moment, {E^) is also shown in Fig. ^. The 
expressions for the second moments converge more slowly in 1/Mb than those for the first 



moments, and the theoretical advice [0 is not to put much trust in the bands they define. 
Consequently we have not shown them in Fig. 0. 

The intersection of the two bands from the first moments determines A and Ai. A 
Ax^ = 1 ellipse is shown. The values obtained are 

A = 0.35 ±0.07 ±0.10 GeV , 



Ai = -0.236 ± 0.071 ± 0.078 GeV^ . 

Here, the first error is from the experimental error on the determination of the two moments, 
and the second error from the theoretical expressions. (Using the information from all four 
bands, first and second moments, the results differ little, both as to central values and as to 
errors.) Note that A and Ai are scheme and order dependent. The values obtained above 
are for A and Ai to order 1/M^, order l3oa1, in the MS renormalization scheme. 

Given this determination of A and Ai, we can use them to improve the determination 
of \Vcb\ from the measured B —>■ X^lv semileptonic width. The expression [l^J^ for the 
semileptonic width, to order PqoI and 1/M|, is given in Eq. ^ 

Tsi = ^%£^0.3689[1 - 1.54f - 1.43/3oft - 1-648^(1 - 0.87f ) - 0.946|^ - 3.185^ 



B 



+ 0-02,% - 0.298|^ - 3.28|| ± 10.47|| - 6.153-§^ ± 7.482-^ 

- 7-4^ ± 1.491^ - 10.41^ - 7.482^ ± 0(1/M^)] . (6) 

For the experimental determination of F^/, we use: B{B Xcdv) = (10.39 ±0.46)% [|TU|, 
tb± = (1.548 ± 0.032) ps p|, tbo = (1.653 ± 0.028) ps 0], /+-//00 = 1.04 ± 0.08 0, 
giving F,; = (0.427 ± 0.020) x 10"^° MeV. 

Combining the measured semileptonic width with the theoretical expression for it, and 
using the determination of A and Ai from the first moments, we find 



\Vcb\ = (4.04 ± 0.09 ± 0.05 ± 0.08) x 10"^ , 

where the errors are from experimental determination of Tgi, from experimental determina- 
tion of A and Ai, and from the 1/M| terms and scale uncertainty in a^, in that order. This 
gives a determination of |V"cf,| from inclusive processes, with a precision of ±3.2%. This result 
depends on the assumption of global parton-hadron duality, with its unknown uncertainties. 

Summarizing, we have measured the first and second moments of the hadronic mass- 
squared distribution in the B meson semileptonic decay to charm, B —>■ Xciiy. We find 
(M|: - M|,) = 0.251 ± 0.023 ± 0.062 GeV^ {{M]^ - M^f) = 0.639 ± 0.056 ± 0.178 GeV^ 
and {{Ml - {Ml)f) = 0.576 ± 0.048 ± 0.163 GeV^ The measurement of (M^ - ) and 
the HQET expression for this moment are used, in conjunction with similar information on 
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FIG. 2. Bands in A - Ai space defined by {M^ — M^) , the measured first moment of hadronic 
mass-squared, and (E^), the first moment of the photon energy spectrum in & — > 57 |Q. The 
inner bands indicate the error bands from the measurements. The light grey extensions include the 
errors from theory. All bands are derived from 0{1/M^) 0{Poa'^) HQET expressions, using the 
MS renormalization scheme. 

the first moment of the photon energy spectrum in 6 ^ 57, to determine Ai and A. These 
in turn are used, along with the B meson semileptonic width, to obtain Vcb- 
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